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We have developed a model for the N N — > NNtttt reaction and evaluated cross 
sections for the different charged channels. The low energy part of those channels 
where the pions can be in an isospin zero state is dominated by A r *(1440) exci- 
tation, driven by an isoscalar source followed by the decay N* — > N(mr)'g^^ llave . 
At higher energies, and in channels where the pions are not in T=0, A excitation 
mechanisms become relevant. 



Pion production in iViV collisions is one of the sources of information on 
nucleon-nucleon interaction and resonance properties. Nowadays, two pion 
production in pp collisions is a subject of experimental research at the CEL- 
SIUS storage ring by the WASA/PROMICE collaboration §. A direct com- 
parison of our theoretical results with the measured total cross sections and 
invariant-mass distributions will certainly provide useful information about the 
mechanisms governing this reaction. 

Our model can also be relevant to understand certain features observed 
in other two pion production reactions; namely, n + p — > d + (irir) and 
p + d — > 3 He + 7r + 7r~ . The first of them was, experimentally studied at 
LAMPF for a neutron beam energy of 800 MeV0, and is also an important 
ingredient of the deuteron spectra measured with a 1160 MeV neutron beam at 
SaturnetJ, where the ABC effect was first observed in the free nuclepn reaction. 
The second is being studied at COSY in the MOMO experiment □. 

Finally, the present model has also repercussions for the reaction p + p — > 
p + p + 7T° since a possibly relevant mechanism can be obtained from p — > 
p + p + TT°n when one of the 7r°'s is emitted and the other is absorbed cl 

In order to build the model, we have closely followed the guidelines of a 
previous model for the ttN — » ttttN reactionu. A relevant finding of that work is 
that, even at threshold, the Roper resonance excitation and subsequent decay 
into Af(7T7r)g^° is a very important mechanism. In our case, the situation 
is similar, but more involved because one has to deal with the NN — > NN* 
transition, which is poorly known. Below, we present a brief description of the 
ingredients of the model. Details can be found in Ref. 7. 

Diagrams (1), (2) and (3) of Fig. 1 can be derived from the lowest order 
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Figure 1: Complete set of Feynman diagrams of our model. 



SU(2) chiral Lagrangians containing both pions and nucleons at tree level 

c = c 2 + c[ B) . (1) 

The expressions for these Lagrangians are given, for instance, in Ref 8. With 
these ingredients we generate only the isovector part of the s-wave ttN am- 
plitude, needed for diagram (3). The much smaller isoscalar part appears at 
higher order. We just consider it empirically 

C = -4tt— dtpi) , Ai = 0.0075 . (2) 

In spite of the fact that these amplitudes do not vanish at threshold, their 
contribution to the total cross section is very small, even at threshold, in most 
channels ( see Fig. 2 ). 

In order to obtain the amplitudes for diagrams (9)-(15) we need the fol- 
lowing phenomenological Lagrangians 

Can* = —^Sjid^T^N + h.c. (3) 

Caa* = —ASAi(di(p X )T^ N + h.c. (4) 

where (T^) and 5a (Ta) are the spin (isospin) 1/2 — > 3/2 transition operator 
and spin (isospin) 3/2 operator respectively,-. The coupling constant /*= 2.13 
is obtained from the A — > Nir partial width la and /a =4/5 Jnntt comes from 
SU(6) quark model. Our T = 1 exchange potential includes ir and p exchange 
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and short range correlations L3. All these terms are quite small except diagram 
(12), which becomes relevant at energies T p > 1 GeV. 

We us^the information about the AT* (1440) properties contained in the 
PDG booku in order to extract the couplings required for diagrams (4)- (8): 

£n*Ntt = ^ N ,a t {d t (}))Ttl> N + h.c. (5) 

£^ A7r = m^lL^ sUd^T^, + h.c. (6) 
ma 

Here, the couplings / = 0.477 and givAn — 2.07 are obtained from the corre- 
sponding partial widths assuming branching ratios of 65 % and 25 % respec- 
tively, and a total width of 350 MeV. 

The Lagrangian for the N* — > TV(7T7r)g=£, al , e decay channel is given in Ref. 
12. After expanding in pion fields one gets 



jCn'Ntttt = -c* ir ^-ipN*<f> 2 4>N - C2j^j^{d f2 d„'ilJN*)(Td tl (f))(Td lJ ct))-if; N + h.c. 

(7) 

In this case, the free parameters c\ and c\ can not be both obtained from the 
partial decay width. They can just be constrained to an ellipse til. In order to 
further constrain these parameters we use the model of Ref. 6 for the reaction 
ir~p — > n + ir~n. The best agreement with the experiment is obtained with 
c\ = — 7.27 GeV^ 1 and c\ = ( Set I ) but the experimental errors are also 
compatible with c* = -12.7 GeV- 1 and c* = 1.98 GeV' 1 ( Set II). 

In this case, the T=l exchange is similar to the one for the A terms, but in 
addition we must consider an exchange in_the T=0 channel. In a recent analysis 
of the (a, a') reaction on a proton target Ej, the strength of the isoscalar NN — » 
NN* transition was extracted by parameterizing the transition amplitude in 
terms of an effective "ct" , which couples to TV TV as the Bonn model a E3 

Conn = g^NNip^, 91nn/4k = 5.69 (8) 
and couples to NN* 

C<tNN* = gcrNN'^N'tfrP + h.C. (9) 

with a strength provided by a best fit to the data: g^jvjv* 1^ ~ 1-33, and such 
that g a NN and g<,NN* have the same sign. 

The amplitudes corresponding to diagrams (4), (5) give a non- vanishing 
contribution at threshold, which is by far the dominant one below 1 GeV for 
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those channels where the pions can be in T=0. The contribution of the T=0 
exchange is an order of magnitude larger than the T=l exchange, which is 
diminished due to the short range correlations. 




Figure 2: Total cross sections for two of the channels, as a function of the incoming proton 
kinetic energy in lab. frame. Solid line, total ( line labelled 1 for set I and 2 for set II ); 
long-short dashed line, N* — » ^{ 7T7T )'g^% ave ) long-dashed line, N* — > An; dash-dotted line, 
A excitation mechanisms; short-dashed line, non-resonant terms from diagrams (l)-(3). The 
partial contributions are calculated with set I. Experimental data are taken from Ref. 14. 



We implement the final state interaction (FSI) in a simplified version of 
the model, meant to work at T p < 900 MeV and for the pp — > ppir + ir~ , where 
we can assume that the total cross section is clearly dominated by the N* 
excitation driven by an isoscalar source. Since the energy of the incoming 
nucleons must be high enough to produce two pions, we neglect initial state 
interaction. For the configuration of the figure, the amplitude is 



M' si = I j0y>M free (q,p N *)MP) (10) 



where 



, P4-P3 „ . P5+P6+P2-P1 . . 

fc = o > P = <1^ o » PN» ~ P5 + PS + P3 , (11) 



For the wave function in momentum space (f>k{P) we use an analytic expres- 
sion obtainjed for a separable non-local potential that takes into account the 
repulsionlia. The parameters are chosen t3 to fit the 1 S'o np phase shifts. 
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Figure 3: Dominant mechanism for pp —* ppir + n~ at low energies with FSI included and 
total cross sections obtained with ( dotted lines ) and without ( solid lines ) FSI for both 
sets I and II. 
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